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Abstract. We have found new anisotropic vacuum solutions for the scale-invariant 
gravity theories which generalise Einstein's general relativity to a theory derived from 
the Lagrangian R 1+d . These solutions are expanding universes of Kasner form with 
an initial spacetime singularity at t = and exist for —1/2 < 5 < 1/4 but possess 
different Kasner index relations to the classic Kasner solution of general relativity to 
which they reduce when (5 = 0. These solutions are unperturbed by the introduction of 
non-comoving perfect-fluid matter motions if p < p on approach to the singularity and 
should not exhibit an infinite sequence of chaotic Mixmaster oscillations when S > 0. 
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There have been many studies of gravity which consider generalisations of the 
Einstein-Hilbert action of general relativity by the addition of terms of higher-order in 
the curvature to the Lagrangian, see refs in [1, 3]. These are of interest for the evolution 
of cosmological models close to the initial singularity and at late times, for the study 
of quantum gravitational phenomena, inflation, and the observational consequences on 
solar system scales. In particular, it is possible to study the structure of gravity theories 
in which the Lagrangian is quadratic in the curvature or an arbitrary analytic function 
of the curvature [4, 1, 5]. Unfortunately, such studies are complexified by the lack of 
exact solutions of the field equations for these theories other than the special cases 
provided by the solutions of general relativity. The field equations are typically fourth 
order and difficult to deal with except for spacetimes with isotropy and homogeneity 
or static spherical symmetry. We have recently considered in detail a different type 
of generalisation of Einstein's theory which derives from a power-law Lagrangian of 
the form R 1+5 , where 5 is constant. This theory reduces to general relativity when 
5^0, and both the flat Friedmann universes [3, 6, 7] and static spherically symmetric 
solutions [6] can be found, tested for stability, and confronted with cosmological and 
perihelion-precession observations to produce tight observational limits on S of [6] 



In this letter we will show that it is also possible to find exact anisotropic and 
spatially inhomogeneous cosmological solutions of this theory of Kasner type. These 
seem to be the first exact anisotropic solutions of higher-order Lagrangian theories of 
gravity and provide a testing ground for a range of interesting investigations. Kasner 
solutions [8, 9] form the building blocks for an understanding of the behaviour of the 
most general known solutions of the Einstein equations and allow us to understand the 
conditions under which chaotic behaviour arises and persists on approach to the initial 
cosmological singularity [10, 11, 13]. They also provide a simple environment in which 
to study quantum gravitational effects like particle production [14], and the conditions 
under which shear anisotropy dominates the expansion of the universe at early and late 
times. The simplicity of these new solutions is appealing and allows for a complete 
understanding of their behaviour as 5 varies. They also provide an instructive context 
in which to evaluate the nature of the general relativistic evolution and to determine 
whether it is typical or atypical within this wide class of theories. 

We consider here a gravitational theory derived from the Lagrangian density 



where 5 and x are constants. The limit 5^0 gives us the familiar Einstein-Hilbert 
Lagrangian of general relativity and we are interested in the observational consequences 
of | <5 1 > 0. There is a conformal equivalence between this theory and general relativity 
with a scalar field possessing an exponential self-interaction potential [15, 16]. 

We denote the matter action as S m and ignore the boundary term. Extremizing 



< 5 < 7.2 x 10 
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with respect to the metric g a b then gives [17] 

5(1 - 5 2 )R 5 ?^ - 5(1 + 5)R 5 ^ + (1 + 5)R s R ab (2) 
R z R 

- l -g ab RR & - g ab 5(l - 5 2 )R 5 ^^- + 5(1 + 5)g ab R 5 ^ = |r a6 , (3) 

where T ab is the energy-momentum tensor of the matter, and is defined in the usual 
way. We take the quantity R 5 to be the positive real root of R throughout. 
We seek solutions of these equations for the Bianchi type I metric 

ds 2 = -dt 2 + t 2pi dx 2 + t 2p2 dy 2 + t 2p3 dz 2 (4) 

where the Kasner indices, Pi, are constants. This metric is a solution of the vacuum 
field equations of the R 1+s theory if and only if the Kasner indices satisfy the algebraic 
relations 



(1 - 5)(p\ + p\ + pi) + 35(1 + 25) = (1 + 25 2 ) ( Pl +p 2 + p 3 ), (5) 
(1 - 5)(p lP2 + pm + P2P3) + 35 2 (1 + 25) = 5(2 + 5)( Pl + p 2 + p 3 ) (6) 

These constraints can be solved to yield two classes of solutions. The first has 
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These are only solved by the isotropic solution with 

, (1 + 25) 

Pl=P2 = P3 = 5 JT^5)- 

This is the zero-curvature vacuum Friedmann universe found by Bleyer and Schmidt 
[2, 3]. The second class of solutions to (5)-(6) is new and generalises the anisotropic 
Kasner universe of general relativity to 



3 

5>=l + 25 (7) 

i=i 

EpH1-45 2 . (8) 

i=i 

Note that when 5 = this reduces to the standard Kasner solution of general relativity 
[8]. The first constraint (7) means that we must have 5 > —1/2 for an expanding 
universe. The second constraint (8) requires —1/2 < 5 < 1/2 for consistency. When 5 
falls outside this range then anisotropic solutions with this simple power-law form do 
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Figure 1. The intervals in which the Kasner indices pi lie can be read off this graph. 
For any value of 8 in the range — \ < S < j a horizontal line is drawn; the boundaries of 
the intervals in which the pi lie are then given by the four points at which the horizontal 
line crosses the two closed curves defined by (9)-(ll). For 5 — these boundaries can 
be seen to be — p ; | and 1, as expected for the Kasner solution in general relativity. 



not exist. These constraints restrict the range of variation of the individual pi as follows 
(assuming with out loss of generality that they are ordered p 1 < p 2 < p^): 

1 + 25-2^/(1 + 25X1-45) ^ ^ 1 + 25 - y/(l + 25)(1 - 45) 

s pi s (yj 

S p2 < l+M + V(l + M)(l-4f) (10) 
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1 + 25- v/(l + 25)(l 
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1 + 25+ v /(l + 25)(l 


-45) 



l + 25 + 2v/(l + 25)(l-45) 

3 -" 3 - 3 • 

These ranges are plotted in Figure 1 for different values of 5. They impose the further 
constraint —1/2 < 5 < 1/4 on the range of allowed values of 5 for which real- valued 
solutions exist. A horizontal line of constant 5 intersects the two closed curves at four 
points which give three line segments covering the allowed ranges of the Kasner indices 
for all — 1/2 < 5 < 1/4. For 5 = 0, the intersects give the ranges for the Kasner solution 
in general relativity: — 1/3 < p\ <0<p2<2/3< 1. 

We note the following features of these solutions, they have a Weyl curvature 
singularity at t = so long as 5 > —1/2. Solutions with 5 < —1/2 are contracting 
universes and have a 'big rip' singularity as t — > oo. When 5 > it is possible for all 
the Kasner indices to take positive values for every positive value of 5. When 5 < we 
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see that p 1 < always and p 2 can also take negative values but the universe expands in 
volume. The measure of the positive values for pi increases as S increases to its maximum 
allowed value. This situation is of particular interest in connection with the behaviour 
of more general Bianchi type IX Mixmaster cosmologies in these gravity theories. In 
general relativity we know that these cosmologies exhibit chaotic oscillatory behaviour 
on approach to the initial Weyl curvature singularity at t = in vacuum [10, 11, 12, 13]. 
In order for an infinite number of chaotic oscillations to occur we need one of the Kasner 
indices to take only negative values. This is the case in vacuum in general relativity. 
When a scalar field is added the Kasner relations change and it is possible for all indices 
to be positive and the oscillatory sequence of permutations of the scale factors will cease 
after a finite number of oscillations [18]. The subsequent evolution will be in general 
of Kasner type. We expect a similar sequence of events in the Mixmaster cosmologies 
in R 1+s gravity theories when 5 > 0. Oscillations of the scale factors as t — > will 
eventually permute the Kasner indices into one of the combinations in which they are 
all positive and oscillations will cease leaving the evolution to proceed in Kasner form. 
By contrast, we expect chaotic oscillations to persist when — 1/2 < 5 < although their 
detailed structure [19] will differ from those that occur in general relativity and will be 
investigated elsewhere. 

The new form of the Kasner index constraints found here differs from the form found 
in general relativity when a massless scalar field is added [20, 23] because a change is 
made to the sum of the indices (7) as well as to the sum of their squares, (8). 

We can also investigate the extent to which this vacuum metric (4) is a stable 
asymptotic solution of the field equations in the presence of perfect-fluid matter as 
t — > as it is in general relativity. We need to determine whether the form of the 
exact solution leads to fluid stresses which grow faster than the vacuum anisotropy 
terms (which are all 0(t~ 2 ( 1+<5 )) as t — > 0). Similar investigations can be made mutatis 
mutandis for the case of fluids with anisotropic pressures [21]. 

Consider fluid motions in the background Kasner metric and assume that the 
material content of the universe is a perfect fluid with equation of state p = (7 — l)p 
with 1 < 7 < 2. The continuity equation is [22, 20, 24] 



where u 1 is the normalised 4- velocity {u a u a = 1), and the momentum conservation 
equation is 



Neglecting the space derivatives with respect to the time derivatives, so that we 
confine attention to scales larger than the particle horizon in the velocity-dominated 
approximation, we have 



d_ 




(12) 




(13) 



= const ; 



,(t-i)/7 = const 
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For relativistic motions, keeping the dominant velocity component (u 3 = u%t 2p:i and 
Uq ~ u 3 u 3 ~ (u 3 ) 2 t~ 2pa ) as t -> gives to leading order p ~ f-7(pi+i*)/(2-7) an d 
■u a ~ ^(pi+P2)(7-i)/(2-7) i Using these asymptotic forms we can now check that all the 
components of which they induce on the right-hand side of the field equations (3) 
diverge more slowly than the vacuum terms, t~ 2 ( 1+s \ as t — > 0. This is the condition for 
the Kasner solution to be unperturbed by the metric effects of the fluid motions. We 
have to leading order that 

T ° ~ pul ~ t- 1 - 2 *-** ; Tl ~ p ~ r ^(i+2*-Ps) (14) 

T| ~ pM2M 2 ~ -d+ 25 -f3) . T 3 „ pU3U 3 „ t -i-2s- P3 (15) 

In the simple case of 7 = 4/3 where the 4- velocity of the fluid is comoving (u, = 5°) 
we require only that p ~ £-2-4<5+2p 3 ,-[i ver g es glower than £~ 2 ( 1+<5 ) and, this requires only 
that ps > 5, which is always satisfied. For general 7 the worst divergence created by 
relativistic motions is in pu s u 3 ~ £-i-2<s-p3 anc [ this is always slower than t~ 2 ( 1+<5 ) since 
P3 < 1. But for stiffer equations of state we have to careful that the assumption that 
the velocities are increasingly relativistic (u° » 1) as t — > continues to hold. This 
requires 7 — 1+25 < p 3 . If this inequality fails then u a u a — > and -u ~ 1 as t — > and we 
need to re-solve (12)-(13) in the approximation where the velocities u a u a go to zero and 
u° ~ 1 because the ultra-stiff fluid makes the motions grind to a halt [25]. In this case 

p ~ £-7(1+2,5) and ^ _ i ( 7 -l)(l+2i)_ g ince we haye _ (M 3 ) 2 t- 2ps ~ t 2( 7 -l)(l+2«5)~2 P3 

this behaviour occurs when 

( 7 -l)(l + 2(J) >p 3 - (16) 

Since u a u a — > as i — > we will always have p » pu a u^ as £ — > and we need 
only check that p diverges more slowly than the vacuum terms 0(t~ 2( - 1+s ^) in order to 
check that the Kasner solution for the metric is unperturbed by the fluid. This requires 
(7 — 2) (25 + 1) < and since 7 < 2 we require 25 + 1 > for the vacuum term 
to dominate. This is just the condition for the universe volume to be expanding as t 
increases, (7), and always holds. 

Hence, the Kasner solutions we have found will provide a good description of 
the general spatially homogeneous perfect-fluid solutions of the R 1+s gravity theories 
in the vicinity of an initial cosmological solution. They may also provide a useful 
approximation to the time dependence of a general inhomogeneous cosmological solution 
to these theories. 

In summary, we have found new anisotropic vacuum solutions for the scale-invariant 
gravity theories which generalise Einstein's general relativity to a Lagrangian R 1+s . 
These solutions are expanding universes of Kasner form with an initial spacetime 
singularity at t — and exist if— 1/2 < 5 < 1/4 but have different Kasner- index 
relations to the classic Kasner solution of general relativity when 5^0. These solutions 
are unperturbed by the introduction of non-comoving fluid matter motions if p < p on 
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approach to the singularity for this range of S, and do not exhibit an infinite sequence of 
chaotic Mixmaster oscillations when 5 > 0. They should provide a simple new testing 
ground for quantum cosmological processes and late-time behaviour in theories 'close' 
to Einstein's general relativity. 
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